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Abstract 

The processes of p-wave Quarkonium exclusive decays to two mesons are investigated, in 
which the final state vector mesons with various polarizations are considered separately. In 
the calculation, the initial heavy quarkonia are treated in the framework of non-relativistic 
quantum chromo dynamics, whereas for light mesons, the light cone distribution amplitudes 
up to twist-3 are employed. It turns out that the higher twist contribution is significant 
and provides a possible explanation for the observation of the hadron helicity selection rule 
violated processes Xci —> 4>4>, ww by the BESIII collaboration in recently. We also evaluate the 
Xbi — > J/iftJ/ip process and find that its branching ratio is big enough to be measured at the 
B-factories. 
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I. INTRODUCTION 



As one of the most interesting fields of high energy physics, the study of heavy 
quarkonium plays an important role in understanding the configuration of hadrons and 
the non-perturbative behavior of QCD. Within the quarkonium physics, recently the 
exclusive double charmonium production gains much attention in both theory and ex- 
periment. The double charmonium production at the LHC experiment [l| provides 
an opportunity to test the prevailing effective theory of quarkonium physics, the non- 
relativistic chromodynamics(NRQCD) |2|, in high energy, and the double charmonium 
production at B-factories J3 1 poses a challenge to the leading order(LO) perturbative 
QCD(pQCD) calculation |^-|6j, which stimulated a series of followup theoretical studies 
7H1I|. 

For the double quarkonium production, since generally the color-octet mechanism 
[2I deduced from the NRQCD formalism is not that crucial as in some other situations, 
for instance in the inclusive quarkonium hadroproduction, it is tempting to apply the 
light cone formalism to the study 9]. At the B-factory, although the center-of-mass 
energy is only about 10 GeV, not very high, it is still much larger than the charm quark 
mass. And hence, the charm quarks in the double charmonium production may move 
on the light cone before hadronization. This fact inspires one to investigate the double 
charmonium production processes at the B-factory energy in the light cone formalism, 
though the form of charmonium light cone wave functions are far from mature. 

At the B-factory, there are a huge number of bottomnium states, and the bottomo- 
nium decays to charmonium pairs are now under investigation in experiment 12] . The- 



oretically, several calculations on this issue are performed in the framework o 
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16], but only a few are carried on in light cone formalism. In Refs. 17H19] the light 



NRQCD 



cone formalism was applied to calculate the process of bottomonium decays to double 
charmonium. The authors of [l^, found that higher twist terms were important and 
the results of Br(j]b — » J /if) + J/if>) are much larger than the leading order NRQCD [jjj] 
prediction, which asks for further studies on quarkonium production and decays in light 
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cone framework. 

In the literature, the p-wave charmonia, the Xcj, exclusive decays to light meson pairs 
have been investigated extensively, by virtue of NRQCD Il3|, the light cone formalism 
with perturbative QCD(pQCD) factorization approach [20], and also the meson loop 



technique 211], respectively. However, the corresponding study in bottomonium sector, 
i.e., the study on the p-wave bottomonia, the XbJ, exclusive decays to light meson pairs 
and double charmonia are very limited. In Refs. [jij, , p-wave bottomonia decays to 
double J/ip are evaluated in light cone formalism at leading order in twist expansion 
and in the framework of NRQCD, but the results from different approaches do not agree 
with each other. 

In this work, we study various processes of p-wave heavy quarkonium exclusive decays 
to meson pairs (VV or PP), in which the light cone distribution amplitudes up to twist-3 
are applied to the final states. Notice that neither relativistic correction in NRQCD, nor 
the twist- 2 contribution in light cone formalism can explain well the BESIII experimental 
data of Xci 'P'Pi UUJ processes, which implies that the helicity selection rule is seriously 
violated. We tend to think that the twist-3 contribution is not merely a higher order 
correction, but rather significant in the description of these processes. 

The rest of this paper is organized as follows. In Section II we present the strategy 
and formalism of the study; in Section III, we perform the numerical calculation and give 
some discussion on the results; Section IV is devoted to a summary and conclusions. For 
the sake of the readers convenience, some of the formulas used are given in the Appendix. 



II. CALCULATION SCHEME DESCRIPTION AND FORMALISM 
A. The initial heavy quarkonium 

To begin the calculation, we first determine the character of initial heavy quarkonium 
via partons scattering process QQ —¥ qq + qq, where Q stands for the heavy quark b 



3 



FIG. 1: Leading-order QCD Feynman diagrams for \QJ ~^ VV or PP process. 



or c, and q for the constituent quarks of final state mesons. The schematic Feynman 
diagrams at the leading order in perturbative QCD are shown in Figure [TJ For initial 
heavy quarks Q and Q, their momenta are assigned as 

P P 

respectively. Here, the relative momentum / between heavy quarks in the center-of- 
mass system(CMS) is much smaller than the total momentum P of the quark-antiquark 
system in the laboratory system(LS). To construct the spin-triplet states xqj, one needs 
to project the colors and spinors of Q and Q to the proper quantum number of the states 



through projection operator by the standard approach in NRQCD 



22|, that is 



3 " 1V2E,(E, + rn Q ) " s/K 

where Ef = P 2 /4 = rriq — I 2 , N c = 3, and l c is the unit color matrix. The spin 
polarization vector satisfies P ■ e = and e ■ e* = —1. 

Generally, in the spin-triplet case the expansion of the matrix element in powers of 



I has the form J4] 

M[QQ(s = 1)] = (A p + B pa r + c paT rr + . . .) 4 ■ ( 3 ) 

Then the matrix element for the spin-triplet P-wave heavy quarkonium states can be 



4 



expressed as: 

M\xqo] 
M[xqi(X)\ 
M[x Q 2(X)} 



(Oi) 



XQO 



1/2 



2N c m Q 

2N c m Q 

(Oi) 
2N c m Q 



B 



pa 



2m Q V2 



e^P A e M (A) 



1/2 



pa 



2 3 



(4a) 
(4b) 

e^(A) . (4c) 



Here, (Oi) XQJ are the NRQCD matrix elements; the tensor J**" reads 



pp pi/ 

4m| ' 



(5) 



and the sums of various polarizations of spin-1 polarization vector and spin-2 polarization 
tensors give 



jy(A)e"(A) = J"", 

A 

2 3 



(6a) 
(6b) 



B. The light-cone distribution amplitudes 



To calculate the hadron's distribution amplitude, the bi-spinors vu of the quark- 
antiquark pairs of final state mesons should be replaced by corresponding light-cone 
)rojectors. The light-cone projectors for light mesons can be readily obtained from Ref. 



231 ]. which is slightly different from what we used in the following by a factor of (— i) in 



convention. The two-particle light-cone projection operator for a light pseudoscalar in 
momentum space up to twist-3 reads 



-h 



< 75 (p{u) - /ip7 5 ( (p p \u) - ia ^ nin 
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+ ivpup 



//■< 



6 dk ±u 



(7) 

Here, p' stands for the momentum of meson; the parameters fip are m^j(m u + md), 
m 2 K /(m u + m s ) and m^ c /2m* in the MS scheme for pion, kion and rj c , respectively; 



</> p (u) and CT are twist-3 distribution amplitudes; n_ and n + are light-like vectors sat- 
isfying n 2 = 0, n\ = 0, n_ ■ n + = 2; the component quarks' momenta are assigned 

as 

fcf = uM + ifei + i;< ) k% = uEn»-k^ + ^-n». (8) 

Note that the term which involves the transverse momentum derivative acts on the hard 
scatter amplitude before collinear limit k\ — up' — uEn- is taken. The asymptotic limit 
of the leading twist distribution amplitude takes the form <p{u) = 6uu. The two twist-3 
distribution amplitudes (fi p and CT can be obtained by solving the equations of motion 
and are determined as <f> p (u) = 1 and 4> a {u) = 6uu. 

In the momentum space, the transverse and longitudinal light-cone projections of 
vector meson reads 

l/r v fv m v (e-n + ) f^m v m v (e ■ n+) \ i ^ {t) 
» = T 2E E ^~ M u ) + — 2E I ~2 <V + II ^ 



W,^ _ ^ d , ^ fv(m 1 + m 2 )^ s Hu) 



zE / dv (Mv) ~ * ») a^nt ^— + (1 - " ) , 



(9) 

k=up' 



and 



M± = -rfvEf±t-Mu) + lfvmv 



? ± g^\u) - Ei_ £dv fa(v) - <#>(„)) 



(10) 

k=up' 



4 <9A; ±(T 

respectively. Here, the transverse polarization vector 

e£ = e" — <■ ( U ) 

In the CMS of the initial quarkonium state, the light-cone vectors of final state mesons 
can be chosen as 

n_ = (1, sin^, 0, cos#) , n + = (1, — sin6>, 0, — cos#) . (12) 
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For the vector meson moving along n_, its polarization vectors are 

cos 9 i sin 9 , 



e + 



6(0) 



(0, 



(o. 



cos # z sin # 

"71" '"71 ' _ 71" ) ' 



^E 2 - m 2 v E sin 9 E cos # , 



my 



(13) 



The higher twist LCDAs are related to the leading twist ones through Wandzura- 
Wilczek relations 



24 1 



ffi o) («) = 2 



L </o 

u 



u 


























If 






"9 




>x(i>) 



0|| (u) 



dv 



dv + u j ^ dr 



dv 



dv 



dv + u I < ^ >± dr 



(14) 
(15) 
(16) 
(17) 



Generally, while the above distribution amplitudes are convoluted with the hard 
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271. In this 



part of a specific process, the troublesome endpoint divergences appears 
work the logarithmic and linear infrared divergences, appearing when the convolution 
integrations are performed, are attributed to certain complex quantities similar as the 



prescription used in Refs. 



u 



25 



/ 

Jo 



-|27J, as 
Inu 



1 



du—^--(X a y 
u 2 



u- 



(18) 



This measure is different from what in the pQCD approach taken in Ref. |20j, where the 
Sudakov form factor is introduced in the convolution integration. Here, the first and 
the second are just the same as given in Ref. {26], and the third one should be less 
than 1 in magnitude, which is smaller than what in Ref. 26] because here the integral 



does not divergent at u = 1. We take the same assumption about X a and Xi as in Ref. 



261 ]. i.e., they are universal to all final states with magnitudes around \n(m,Q/ Aqcd) and 
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tjiq/Aqcd, respectively. And, hence they can be constructed as 

X a = (l + p a e^)M-^), X I= (l + pie <*)J22-. (19) 

The values of p and (f) are such constrained that p < 1 and is arbitrary In this 
work we take p a = Pi = 0.5 and (fi a = <fii = 90° by default, which is slightly different 
from what used in Refs. j^, Q in B and B s decays. The asymptotic amplitudes 
4>{u) = (j>a(u) = <p±{u) = (j)\\(u) = Quu are taken in the calculation. To test the parameter 
dependence, we let pi and (pi to vary the same way as p a and 4> a , respectively. Note that 
all these procedures in dealing with the divergences lead the calculation results to be 
model dependant. 



C. The decay width 

Once one knows the distribution amplitude, the decay widths can be readily obtained. 
The general polarized decay width reads: 



r J - 1 1 [ dcos6 \M J I 2 (20) 

iAl ' A2 - 2M XQJ 87r2J+l\/ 1 Ml n \_ x 2 \ M ^~ ^ 



1 1 / 4m 2 f 1 d cos 6* , 

IXQJ 87r2j+iY ~Iq^ 

Here, m is the mass of light meson and J stands for the quantum number of total 
angular momentum. In case the final states are identical mesons, a statistical factor of 
1/2 should be added to the width (f20j) . While final states being vector mesons, the total 
width then is: 

Ktai = r o J ,o + 2lt + + 2Tj t _ + 4r J >+ . (21) 



D. The calculation procedure 

In our calculation, the computer algebra system MATHEMATICA is employed 
with the help of the packages FEYNCALC [28|. FEYNCALC is used to trace the 
Dirac matrices. The derivatives respect to the momentum are performed by means of the 
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MATHEMATICA code developed by ourself, and it has been tested before applying 
to the real calculation by recalculating the photonic decay rates of the Xco and x& at the 
leading order and getting an agreement with Ref. Jj]. The final analytical results are a 
little bit too lengthy to be listed here. The analytical expressions for B pa in equation (4) 
for various processes with different final states (VV or PP) are presented in the appendix 
for the sake of the readers convenience while they use. 



III. NUMERICAL CALCULATION AND DISCUSSION 
A. Input parameters 

Before carrying out numerical calculations, the input parameters need to be fixed. 
As for the NRQCD matrix elements, we take those used in Refs. 29h32| . 



-\( Xc0 \^(--T>.a) X 



1.1 



1 , / , , + / i it 
-\{XmW(~-V-a)x 



The leading order running coupling constant 



o)l 2 

0)| 2 

o)l 2 

o)l 2 
o)l 2 
o)l 2 

47T 



0.051 GeV 5 
0.060 GeV 5 
0.068 GeV 5 

2.03 GeV 5 , 
2.03 GeV 5 , 
2.03 GeV 5 . 



(22) 



6 ln(/i 2 /A^ CD ) 

is employed with 6 = 33 ~ 2 "' / . In numerical evaluation, we take the quark flavor number 
rif = 4 and Aqcd = 0.225GeV; the interaction scale n = 2rrib for the bottomonium 
decay processes, and /x = 2m c for the charmonium decay processes; the up and down 
quarks are taken to be massless, while the strange quark mass m s = O.lOGeV, the 
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charm quark mass m c = 1.4 ± 0.2GeV in NRQCD calculation, and the bottom quark 



mass mb — 4.8 ± O.lGeV; the masses of quarkonia are obtained from the PDG 



which are M, 



A 60 



9.859GeV, M. 



M Xcl = 3.510GeV, and M Xc2 



Xbl 

3.556GeV . 



9.892GeV, M 



Xl/2 



9.912GeV, M 



3.414GeV, 



The input parameters for final state vector and 



in Table which are obtained from Refs. [18! 
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pseudoscalar mesons are presented 
]. For r] c , the charm quark mass in 



the MS scheme is taken to be m^ IS = 1.2GeV in light cone calculation [19]. 

TABLE I: Input parameters for final state mesons. 





P 


K* 








m v [MeV] 


770 


892 


782 


1020 


3097 


/y[MeV] 


205 ±9 


217 ±5 


195 ±3 


231 ±4 


416 ±5 


/f[MeV] 


160 ± 10 


170 ± 10 


145 ± 10 


200 ± 10 


379 ± 21 




7r + ("7r~) 


K+(K~) 


Vc 




m P [MeV] 


139.6 


493.7 


2980 




/ P [MeV] 


130.4 ±0.2 


156.1 ±0.8 


373 ± 64 




//p[MeV] 


1500 


1700 


3700 





B. Numerical results and Discussions 

By virtue of the formulae provided in section II and combined with tensors given 
in appendix, one can readily get the analytical decay widths for final state mesons 
with specified helicity. After substituting the input parameters in preceding section to 
the analytical expressions, the numerical results are obtained. The magnitudes of decay 
widths for various processes, including helicity decay widths for final state vector mesons, 
are presented in Tables II to IV. 

In our calculation, the uncertainties are estimated as follows. The first one comes 
from the uncertainty of heavy quark mass tuq. The uncertainty in a s (2mg) induced 
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TABLE II: The polarized decay widths for Xbj — > VV in unit of eV. 





Tn nfeVl 
U, u L J 


r+,+[eV] 


- 1 - +, U L 1 - v J 


T + \eV] 

x +, — L c v J 


Xbo 




9 «n +0 ' 48+0 ' 09 

—0.39—0.09 


rt nr ,„±n oo^-i-o 07Q 
u.uuy_ .002-0.009 








Xw - 




c 9n +0.96+0.18 
J - zu -0.78-0.18 


n n-i n + 006+0 158 

U.U±O_ .004-0.018 








Xbo ~ 




fi otr+0.12+0.38 
u.o<J_0.io— 0.14 


n no ^ +0 01 0+0 31 5 

u - uo ^-0.008-0.034 








Xbo -> 




fi or+0.12+0.38 
u.o<J_0.io— 0.14 


n no ^+0 010+0 315 

u " uo ^-0.008-0.034 








Xbo 


— > WW 


9 i o+0. 39+0.07 
Z-io -0. 32-0.07 


nn __L(] 009-1-0 OfiS 
U.UU/_ .002-0.007 








Xbo 


-»• # 


O QO+0. 71 + 0.41 

°- yo -0.59-0.15 


r\ r\n A +0 01 0+0 306 

u.uo<±_ .007-0.034 








Xbo -> 


T / 1 Til 

J/tpJ/4> 


iu - ' u -l. 10-6.79 


_+3 22+93 96 
lJ--J-O_2.46-10.89 








Xbi 










1 2+°-o + ? 5 o 9 x 1(T 3 

—0.2—1.2 




Xbi - 










9 4+0.6+51.8 1f1 -3 
z -^-0.4-2.4 x iu 




Xbi - 










u -1.1-5.9 




Xbi ->• 


K* + K* 








u - J — 1.1— 5.9 ^ xu 




Xbi 


— > OJUJ 








1 2+°- 3 + 24 ; 9 x 10- 3 

- L --^— 0.2— 1.1 ~ iu 




Xbi 










J —0.6—3.3 




Xbi ->• 


Til Til 

J/i>J/4> 








n r-| +0.11+10.58 
0.09— 0.25 




Xbi 


—ffjp 


n i 4 +0.03+0.03 
u-i -0.02-0.03 


Q cn+0.14+6.80 lf) - 
u - ou -0. 11-0.41 x iu 


-3 


n n9_ +0 - 006+0 - 145 

u,u ^ -0.005-0.023 


n fi1 +0.12+0.17 
u - ux -0.10-0.10 


Xb2 - 




n oo+0. 06+0.06 
U.ZO_ 0.04— 0.06 


n +0.28+13.60 1f1 - 
1 - u -0.22-0.82 x iu 


3 


n n/i«+°- 012 + - 290 
u.u^±o_o .010-0.046 


1 99+0.24+0.34 
l-ZZ_o 20— 0.20 


Xb2 - 




n q R +o.07+o.io 

U - OD -0.05-0.09 


9 94+0.65+30.63 i n 
z - z ^-0.50-1.80 x iu 


-3 


n ns9+°- 019 + - 485 

U.UOZ_ .015-0 .077 


1 R c;+0.33+0.59 
l-DO_o.26-0.34 


Xb2 ->• 


K*+K*- 


n q R +o.07+o.io 

U - OU -0.05-0.09 


9 9/1+0.65+30.63 v i n 
z - z ^-0.50-1.80 x iu 


-3 


n ns9+°- 019 + - 485 
u.uoz_ .015-0.077 


1 fi r+0.33+0.59 
l.UO_o .26-0.34 


Xb2 


— > U10J 


n i O+0.03+0.04 
u - io -0.03-0.03 


u - DD -0. 15-0.59 x iu 


-3 


n nq9+°- 007 + - 188 
U.UOZ_ .006-0.030 


n 7 o+0.15+0.22 
u - '°-0.12-0.13 


X62 


-> # 


n 94 +0.04+0.09 
u - Z4 -0.04-0.09 


9 7 o+0. 80+30.10 1fl - 
z -' 5 -0.61-2.71 x iu 


-3 


u - u ' '-0.014-0.076 


1 r 7 +0.31+0.65 
i - ' -0.25-0.32 


Xb2 ->• 




o oi +0.77+4.84 
°- o± -0.63-3.81 


9 c 9 +0. 73+13.86 
z - JZ -0.56-2.14 




r 4^+1.19+36.39 
d -^°-0.97-4.94 


r-, Q-i +12.46+140.26 
di - yl -9. 85-33.66 



by itlq has also been taken into account in the evaluation. The second one comes from 
the uncertainties of parameters p and 0, and is evaluated by taking the decay width 
as function of p and 0. These two sources of uncertainty are the major ones in the 
calculation of this work. 
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TABLE III: The polarized decay widths for Xcj —> in unit of keV. 





I o, o[keVJ 


I +i + |keVJ 


I +, o[keVJ 


l +, — [KeV J 


XcO 


o n 


o ofi +9.63+0.99 
-J- y 0-2. 62-0.90 


n o-i +1.84+2.81 
U " ii -0.25-0.31 






XcO - 




7 Q9+19.26+1.98 
'• yz -5. 24-1.80 


n fi9+3.68+5.62 
u ' uz -0.50-0.62 






XcO - 


■>■ K K 


o fi7 +19.17+2.89 
°- u( -5.66-2.56 


i 09+7.93+12. 16 
1 - OZ -1.07-1.32 






Xco ->• 


if* + A* 


o fi7 +19. 17+2.89 
°' u/ -5.66-2.56 


1 09+7.93+12. 16 
1 - OZ -1.07-1.32 






XcO 


— )■ WW 


o 1fi +7.58+0.74 
°- iu -2.08-0.68 


n 9fi +1.59+2.71 
u - zu -0.22-0.26 






XcO 


-> 00 


a 91 +7.27+2.36 
^• z± -2.64-0.71 


1 n o+6.56+9.95 
1 - uy -0.89-1.08 






Xcl 


n n 






U.U6t_ Q .028-0.030 




Xcl - 








n n7/i+°- 28 + - 32 
u - u '^-0.056-0.060 




Xcl - 








n nQ7+°- 35 + - 43 

u - uy ' -0.72-0.074 




Xcl ->• 


K* + K* 






n nQ7+°- 35 + - 43 

u - uy( -0.72-0.074 




Xcl 


— )■ WW 






n n9/i+°- 089 + - 102 

U.U_^_ .017-0.020 




Xcl 


J. i 






n -, -, n +0.358+0.476 
U.11U_ .081-0.088 




Xc2 




n Kc+l. 79+0.59 
U - OD -0.39-0.41 


n n9Q+°- 174 + - 271 
U.uzy_ 02 4_o .022 


n 1 1 +0.42+0.65 
u - ii -0.08-0.11 


9 Qo+11.42+2.62 
z.yo_ 2 15 _i 32 


Xc2 - 


-> p+p~ 


i i 9 +3.58+1.18 
i - iz -0.78-0.82 


n n^8+°- 348 + - 542 

U.UdO_ 048 _o 044 


n 99+0.84+1.30 
u - zz -0. 16-0.22 


r o fi +22.84+5.24 
O.OU_ 4 .30-2.64 


Xc2 ~ 


» if* if* 


i ^7+4.43+1.76 
J--^' -1.04-1.24 


n 1 qo+0. 833+1. 109 
u - ±oo -0. 112-0.108 


n q/,+1.20+2.00 
u -°^-0.25-0.25 


QC -+36.90+10.21 
- J,J -6.66-4.49 


Xc2 -> 


K*+K*- 


1 , 7 +4.43+1.76 
1.04— 1.24 


n 1 qo+0.833+1.109 
U.100_g 112-O.IO8 


n q^+1.20+2.00 
u - J ^-0.25-0.25 


QC -+36.90+10.21 
°- yo -6.66-4.49 


Xc2 


— ► WW 


n 7 9 +2.27+0.77 
u -' z -0.51-0.53 


n n-?Q+°- 232 + - 360 
u.u«jy_ .031-0 .022 


n 1 .+0.53+0.83 

"•-■-^-o.io-o.io 


q 70 + 14.8 + 3.45 

O. 'o_2.78-1.74 


Xc2 


-+# 


1 1 r+3.10+1.70 

^^-o.si-i.io 


n 1 cr+0.938+1.090 
U.ldO_ 126-0 .147 


n 9C-+0.80+1.55 
u - zd -0. 18-0.24 


r9+35.46+10.58 
°' dz -6.35-4.21 



The polarized decay widths and r + come solely from the twist-3 contribution, 



which therefore do not exist in Ref. 



181 ] where only the leading twist distribution is 



considered. The twist-3 is the leading contribution for xqi ~ ► W processes. In the 
literature, the polarized decay width r + of xqi from the contribution of QED correction 
had been taken into account [13], however we find the QED contribution is several 
orders less than what from QCD, and hence are negligible. From Tables [Til and II III we 
notice that the twist-3 contributions T + + and r + are prominent for polarized widths of 
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TABLE IV: The decay widths of XQJ ~~ PP processes. 







VcVc 


K+K~ 


7T+7T 


nXbj)[eV] 


Xbo 

Xb2 


20.901|S1I;^ 

A i 7 +0.74+0.69 
-0.62-0.65 


1 7C-+0.31+0.06 
L - ' °-0.27-0.06 

n oc+0.06+0.01 
u - o,J -0.05-0.01 


n orr+0.16+0.02 
u - oo -0.13-0.02 

n i 7+0.03+0.00 
U- 1 ' -0.02-0.00 


r( XcJ )[keV] 


XcO 
Xc2 




9 nc;+7.59+0.77 
^• yd -1.97-0.70 

n 7fi +1.95+0.20 
u - ' u -0.51-0.18 


i kq+4.01+0.30 
i -°°-l. 03-0.28 

n OQ+1. 03+0.08 
u - oy -0.26-0.07 



Xbo/Xb2 ->■ J/i>J/if> processes. While for processes Xco/Xc2 -> and Xda/Xc2 -> 00, 

though the twist-3 contributions are not predominant, they are still important. 

Notice that the higher twist contributions are suppressed by powers of and 
in dealing with the infrared divergences appearing in the next-to-leading contributions 
in twist expansion, factors of ln( A ^ g ), In or A ^ D are induced via the regu- 

larization procedure ffl8|) . Except for polarized decay width T + , which comes from 
the interference of twist- 2 and -3 terms and hence has the power suppression of ^ at 
the amplitude level, other polarized decay widths from twist-3 contributions all have 
power suppression of (^) 2 at the amplitude level as expected. The r 0> contains only 
m ( A^c D ) f° rm , the r+j o contains both ln( A ^ g ) and In 2 ( A ^ D ) forms, while r+ ( + and 
r +) _ contain all three kinds of regularized forms. So, as m v ~ Aqcd, the convergence 
in twist expansion is hold well by power counting. 

Of the processes Xba/Xb2 — > J /if) J /if), the twist expansion fails. For xw ~ ► J/if>J/ip, 
the polarized decay width r +j + , which comes merely from the twist-3 contribution, is 
even bigger than r , o, mainly the twist- 2 contribution. For Xb2 — > J/tf)J/if) process, the 
twist-3 contribution r + o is also larger than r 0i o- This is not a big surprise, though 
mj/i/j is smaller than m XbJ , it is still much larger than Aqcd in comparison with other 
light mesons, which spoils the twist expansions. Hence, the perturbative calculation on 
Xbj —> J/ipJ/if> in this work can be treated as a qualitative estimation. 

To estimate the branching ratios we use the following expressions for the total bot- 
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tomonium decay widths [l8|, |3_4j : 

3C- ^ 



XbO 



Xi)2 



' F „„ 2 (^P ) 

-7ra 3— 



61 



C F a\ 

N r . 



n f^2 (Ob) 



(Op) 
m\ 

587 

~54~ 



3 



mi 



3 ™ s 



■AO*) 



0.68 MeV, 
0.20 MeV, 



317 
288' 



-7T 



c A 



16 
27 



9 2m b 



(0»> n/ 



(23) 

(24) 
2 (0 8 ) 



0.09 MeV. 



(25) 



The total decay widths for Xcj are obtained from PDG [33J, i.e., T Xc0 = 10.3MeV, r Xcl 



0.86MeV, and T Xc2 = 1.97MeV. 

The branching ratios of different processes are listed in Table V. They are evaluated 
by taking the center value of each unpolarized width and divided by the total decay 
widths of xqj- 

TABLE V: The branching fractions of various processes of p-wave heavy quarkonium exclusive 
decays to double mesons. 



P°P° 



K*K* 



coco 



<j><t> J/^J/i> 



Xbo 
Xbi 



3.9 x 1(T 6 
4.8 x 1(T 8 



9.3 x 1(T 6 

2.4 x 1(T 7 



3.2 x 1(T 6 
4.8 x 1(T 8 



5.9 x 1(T 6 

1.3 x i(r 7 



4.9 x 1(T 5 

2.i x i(r 5 



Xb2 


7.3 x 1(T 6 


2.0 x 10~ 5 


9.3 x 10~ 6 


1.8 x 10~ 5 


6.7 x 10~ 4 


XcO 


4.4 x 10~ 4 


10.9 x 10~ 4 


3.6 x 10~ 4 


6.2 x 10~ 4 




Xcl 


1.7 x 1(T 4 


4.5 x 10~ 4 


1.1 x 10~ 4 


5.1 x 10~ 4 




Xc2 


35 x 10~ 4 


11 x 10~ 3 


45 x 10~ 4 


98 x 10~ 4 






VcVc 


K+K- 


7T + 7T~ 






XbO 


3.1 x 10~ 5 


2.6 x 10~ 6 


1.3 x 10~ 6 






Xb2 


2.i x i(r 5 


1.8 x 10~ 6 


8.5 x 10~ 7 






XcO 




2.9 x 10~ 4 


1.5 x 10~ 4 






Xc2 




3.8 x 10~ 4 


1.9 x 10~ 4 
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For comparison, in Table fVTI we also present the total decay width of Xbj — > J/tftJ/ip 
obtained in other calculations. The second and third columns give what we obtained in 



this work. The fourth column shows the NRQCD results of Ref. [131]. including QED and 
relativistic corrections. The last column presents the calculation results in lieht-cone 

n 

formalism up to twist- 2 [18]. 

TABLE VI: The total width of \bj J/ipJ/if) process in various calculation methods. 

[eV] T[twist-2] T[twist-3] r-NRQCD[13] T-LC twist-2[18] 

Xm^J/^J/^ 21.69±|;|| 33.05^^^09.27 5.54 79. ± 3.1 ± 25. ±31. 

Xbi^J/^J/^ 2.04l° : ^+f Q 32 9.04 xlO- 7 

Xb2 -»■ J/ipJ/fp 25.84^; 7 J 134.51±||;||±|| 8 3 ;f 4 10.6 270. ± 41. ± 93. ± 110. 



From the results in Table IVT1 we notice that the light cone formalism generally yields 
more than what from the NRQCD calculation for J/ip production, and the results ob- 



tained in Ref. 



13| are much smaller than ours, especially for the Xbi J/i^J/i^ process. 



Further experimental measurement on Xbi J/ipJ/ip process may provide test about 
these two mechanisms. On the other hand, within the light cone formalism the results 
obtain in Ref. [3| are much larger than ours. We have compared the analytical result of 
the twist-2 part, find that there is a redundant factor of two in Ref. 18j. Furthermore, 
the different choices of input parameters may also induce certain uncertainties. 

The comparison of our calculation with experimental results is given in Table IVII1 
We can see from the table that our results in general agree with the experiment in all 
Xcj VV processes within an order of magnitude, while the NRQCD calculation in 
Ref. [3] for Xci VV processes gives results several orders smaller than the BESIII 
measurements. Except for Xc2 — > K + K~, the calculated Xco/Xc2 to nn and K + K~ 



decay widths are about an order smaller than the PDG[33j] data. To be noted that of 
all our calculation results, the Xcj — > 4>4> processes possess larger branching ratios than 
that of the Xcj — > processes, while experimental data do not tell so all the way. This 
is an open question leaving for further study, and could be partly understood as the 
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TABLE VII: Branch ratios of Xcj ~^ W an d Xcj ~^ PP-> from BESIII experiment 36], the 
PDG data [3^, and our calculation. 

Br(twist-3) BESIII [36] PDG [33J 

Xco -> # (6-2l4 9 3 ^ 2 78 6 ) x 10~ 4 (8.0 ± 0.3 ± 0.8) x 10~ 4 (9.2 ± 1.9) x 10~ 4 

Xd -> # (54± 46 8 ™) x 10- 4 (4.4 ± 0.3 ± 0.5) x 10~ 4 

Xc2 -> # (98t 4 ^+ 459 ) x 10- 4 (10.7 ± 0.7 ± 1.2) x 10- 4 (14.8 ± 2.8) x 10~ 4 

Xco -> (3-6l2 ii"i 6 2°) x 10~ 4 (9.5 ± 0.3 ± 1.1) x 10~ 4 (22 ± 7) x 10" 4 

Xd -> ™ (l-llasiag) x 10_4 ( 6 -° ± °- 3 ± °- 7 ) x 10 ~ 4 

Xc2 -> (45l^f) x 10" 4 (8.9 ± 0.3 ± 1.1) x 10~ 4 (19 ± 6) x 10" 4 

Xco -»• K*K* (10.9± 34 6 °+ 2 6 - 4 ) x 10- 4 (17 ± 6 ± 1) x 10~ 4 

Xco -»• i^+iT- (2.9lo;rlo:?) x 10" 4 (6.10 ± 0.35) x 10~ 3 

Xco -> TTvr (l-5±f;g±g;l) x IP' 4 — (8.4 ±0.4) x IP' 3 

Xc2 -> (10.6j^g^ 8 6 - 4 ) x 10- 3 (2.5 ± 0.5) x 10- 3 

X c2 -> K+K- (3.8±|:f±5;g) x 10" 4 (10.9 ± 0.8) x 10" 4 

Xc2 -»■ vrvr (2.0if;|tg;f) x 10~ 4 (2.39 ± 0.14) x 10~ 3 

difference in the selection of nonperturbative parameters p and <p in various xqj decay 
processes. 



IV. SUMMARY AND CONCLUSIONS 



In this work, various exclusive processes of p-wave Quarkonium decays to two mesons 
are investigated, where the final state vector mesons in different polarizations are con- 
sidered separately. For light mesons, we expand the light cone distribution amplitudes 
up to twist-3, which induces certain end-point singularities in the calculation. To hurdle 
these singularities, the prescription used in dealing with the B meson decay problems is 
employed, in which two complex parameters are introduced. Though these parameters 
induce some uncertainties in our calculation, we believe our result should be reliable as 
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an order estimation. 

The Xci and xti to VV decay processes are induced entirely by the twist-3 contri- 
bution in distribution amplitude expansion, which are firstly estimated in this work. 
It turns out that the higher twist contribution is significant and provides a possi- 
ble explanation for the observation of hadron helicity selection rule violated processes 
Xci 00 an d oju observed by the BESIII collaboration in recently. We also evaluate the 
Xbi — > J/4 ) J/'*P process and find that its branching ratio is big enough to be measured 
at the LHC and B-factory experiments. 

Note that after this work has been done, the Belle Collaboration releases new mea- 
surement results on the p-wave spin-triplet bottomonium decays to double charmonium 



12| . Their obtained upper limits of the branch ratios Br(xbo — > J/^J/ip) < 7.1 x 10 5 , 
Br(xbi -> WW) < 2.7 x 1(T 5 , and Br( X b2 -> WW) < 4.5 x 10~ 5 at the 90% 
confidence level are compatible with our calculation, that is: Br(xw — > J/4 ) J/'>P) = 
(0.7 ~ 35.9) x 10- 5 , Br(xbi ->■ J/tpJ/4>) = (0.7 ~ 44) x 10~ 5 , and Br( Xb 2 ->■ WW) = 
(6.9 ~ 295) x 10~ 5 , while taking the uncertainties in both theory and experiment into 
account. 
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Appendix A: tensor structure of the amplitude 



For vector mesons in the final state, the tensor B^ v reads: 

ftlXV = T 0( B »V + Qp, + ^ + ^ (Al) 



BT = f v 2m v £l ■ Y 2 ' "~ {B kl gT + B k2 (nV + nV) + B k3 {n»_n\ + <r£) 

m Q 

fT^ 2 

+B M gT} + ^^{Bk^lAn- - n\)e 2 • n_ + e^(n v + - n^e, • n_] 
m Q 

+B k& {nf i _e v Ll e 2 ■ n_ + n%e\_ 2 Ei ■ n + ) + B kl {n%e\ x e 2 ■ n_ + nte^ei ■ n + )} 
+ft 2 {B k& {el 1 e\ 2 + el^) + B k9 e ±1 ■ e ±2 <?f } , (A2) 



+C k3 (n!Ln v _ + n+n+)e_u • £±2 + CkA{n+n v _ + n^n+)e± 1 ■ e± 2 
+C fc5 [e^i(n- - <)^2 • n_ + - • n+] 

+Cfc6( n - £ ±i £ 2 ■ rc- + rc+ £ i.2 £ i • n+) + Cfc7(n+£^ 1 £ 2 • ft- + n^e^si ■ n+) 
+C k8 Ei ■ n + e 2 ■ n-g^ + C k o,E\ ■ n + e 2 ■ n^(n^ 1 _n l '_ + n%n\) 

+C kl0 e 1 ■ n + e 2 ■ n_{n%n v _ + n^n u + )} , (A3) 

f2 2 

D"" = ^^{£> fc i(£^£l 2 + e v ±1 e^) + D k2 e ±1 ■ e ±2 g^ + D k3 (e ±1 ■ e ±2 )g^ 

m Q 

+D k 4S ±1 ■ e ±2 (n%n\ + n^_n v _) + D k5 e ±1 ■ e ±2 (rf + n v _ + n^_n v + )} , (A4) 



g M z/ _ -Jxlx. — — {^(e^e^a + e u ±1 e± 2 ) + S^fe^e^ + £ ±i e ±2) 
m Q 

+S k3 e±! ■ e± 2 g^l + S k4 e ±1 ■ e± 2 (ntn u _ + n%n\ - - n%n v _) 

+Sfc5[£i • - n+)e^ 2 + e 2 ■ n_(n+ - 

+ 1 S fc6 «e i 1 / £2 • n_ + nle^i • n+) + S^nt^ ■ n_ + • n + )} . (A5) 
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Here, T° = 2 ,„ 7^ = = and of" = q M1/ — " "+ — " + " comes from the derivative with 

lQ\<2N£y'N c m- i ri 3 ± 3 1 2 



respect to the transverse momentum, such as = of. 
The coefficients in above expressions are: 

B kl = -l{(2(con 2 - 1)X 2 + 21n(2)(l + con 2 + (1 - con 2 )X t ) - AX a + tt 2 } ,(A6) 

8 

5 fe2 = "^{2(X a 2 + (2 + 4con - 2con 2 )X a - ln(2)X l - (1 + ln(2) 

-4conln(2) + 2con 2 ln(2))) + (con 2 - Icon - \)ix 2 } , (A7) 

1 



B k3 = ™ {72 ln(2)X, - 72X 2 - 144(1 - con) 2 X a + (72 - 288con 
128 



2 — ±^^1 

+144con 2 ) ln(2) - 1 + 36(1 - confix 2 } , (A8) 

B k4 = l{2(con - 1) 2 X 2 - 4(con + 2)X a + 2(1 - con) ln(2)X l 
8 

+6(1 + con) ln(2)} , (A9) 
B k5 =l(l- CO n)(2X 2 -7T 2 ) , (A10) 

B k6 = -^{2(conX 2 - 4(1 + con)X a - 2(1 + con)(X i - 5) ln(2)) 

+ (3con + 4)vr 2 } , (All) 

B k7 = ^{con(-2X 2 - 8X a + 4(X, + 1) ln(2) + tt 2 ) - 8X a 

+4(X z + l)ln(2)} , (A12) 



9tt 2 

B k8 =-— , (A13) 



9tt 2 

B k9 = — , (A14) 
C kl = "^{lOA 2 - 8X a + 2 ln(2)X, + 5tt 2 + 6 ln(2)} , (A15) 
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Ck2 = g {6X 2 - 2X a + vr 2 - 2 ln(2)} , 
: ^{-2X 2 + 2X a -41n(2)Xi + 37r 2 -44 + 21n(2)} , 
-^{10X a 2 - 6X a - 4 ln(2)X, - 44 - 2 ln(2) + 5tt 2 } , 
C k5 = 1{2X 2 + 3tt 2 - 8} , 

C fc6 = -^{2X,ln(2)-4-ln(2)} , 

9X a 91n(2)X, 451n(2) 
C fc7 --- - +^g— , 

_9vr 2 

^/c8 — -g— j 

C fc9 =-^(tt 2 -16), 
C kw =-^(tt 2 + 16), 
£) fc i =^{2X a 2 -2X a - ln(2)^ + 31n(2)}, 

o 

D k2 =^{X,-l}ln(2) , 

A 3 =-^{X 2 -X a + ln(2)}, 

'*4 =-^{X 2 + X a -21n(2)X i -2 + ln(2)} , 
20 



D k5 =-{7X a 2 -5X a -2X i ln(2)-2 + 71n(2)}, (A29) 

S kl = ^{4 ln(2)X t - 2X 2 a ~ ^ a + tt 2 + 4} , (A30) 

S k2 =-^{-8X a + 7r 2 + 2 + 81n(2)} , (A31) 

S k3 =-^{4X a + (X,-5)ln(2)} , (A32) 

S ki =^{6X 2 -16X a + 7r 2 + 161n(2)} , (A33) 

S k5 =-^{7r 2 -2X a 2 }, (A34) 

S k6 = ^{2(A 2 - 4X a - 2 ln(2)X, + 6 ln(2)) + 3tt 2 } , (A35) 

S k7 =-^{4(X z + l)ln(2)-X 2 -8X a + 7r 2 } , (A36) 

with con = 1 - M !? 1+m2) - 

For pseudoscalar mesons in final state: 

W u = f 2 { }4^M kl g^ + M k2 g^ + M k3 (n^+nln^ 
m Q 

+M k5 g^ + M k6 (ntn u _ + «) + M kl {n^n\ + n%n v _)}. (A38) 
Here T° is the same as the vector meson case in above, and 

M kl = 6{ln(2) - X a } , (A39) 

M k2 =2{X a -ln(2)} , (A40) 
21 



M k3 = \{2X 2 a - 8X a + 2X t ln(2) + 1 + 6 ln(2)} , (A41) 

o 

Mm = ^{X 2 + (X,-l)ln(2)}, (A42) 
9tt 2 

M k5 =—, (A43) 
9tt 2 

M fc6 =9-—, (A44) 
M k7 =-A{7r 2 + 16}. (A45) 
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